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The minimal Type I see-saw model cannot explain the observed neutrino masses and the baryon 
asymmetry of the Universe via hierarchical thermal leptogenesis without ceding naturalness. We 
show that this conclusion can be avoided by adding a second Higgs doublet with tan /3 > 4. The 
models considered naturally accommodate a SM-like Higgs boson, and predict TeV-scale scalar 
states and low- to intermediate-scale hierarchical leptogenesis with 10 3 GeV < Mn 1 < 10 s GeV. 
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I. INTRODUCTION 

The discovery of a neutral Higgs boson at the Large 
Hadron Collider (LHC) (TJ [2j has strengthened the case 
for the standard model (SM) paradigm of spontaneous 
electroweak symmetry breaking: a scalar doublet < f > 
gains a vacuum expectation value (vev) (<f>) = v/A ss 
174 GeV by virtue of the potential 

V SM = /r 1 2 $ t $ + A($ t $) 2 , (1) 

with fj, 2 < 0. As well, the measurement of neutrino oscil¬ 
lations EHE] suggests that the SM should be extended to 
incorporate neutrino masses. A straightforward way to 
achieve this is to include three right-handed neutrinos. 
Then gauge invariance allows two extra renormalisable 
terms to be added to the Yukawa Lagrangian, 

-A C Y = + I'MArYvr + h - c -> ( 2 ) 

where II = {vl,£l) T , $ = and Mj are the right- 

handed neutrino masses. The SM extended in this way 
is what we refer to as the minimal Type I see-saw model 

[M2- 

If y v v < Mi then the minimal Type I see-saw provides 
an elegant explanation for the smallness of the neutrino 
masses. After electroweak symmetry breaking, the neu¬ 
trino mass matrix is given by the see-saw formula 

v 2 

m v = Iv^M Vv ! (3) 

where T>m = diag(Mi,M 2 ,M 3 ), suppressed by the pre¬ 
sumably large right-handed neutrino mass scale. 

Minimal Type I see-saw also provides a mechanism to 
reproduce the baryon asymmetry of the Universe (BAU). 
Fukugita-Yanagida hierarchical thermal leptogenesis El 
proceeds via the CP-violating out-of-equilibrium decays 
of the lightest right-handed neutrino N lt creating a lep¬ 
ton asymmetry which is reprocessed into the baryon sec¬ 
tor by the electroweak sphalerons. Successful hierarchi¬ 
cal thermal leptogenesis is possible when the Davidson- 
Ibarra bound (ensuring enough CP asymmetry in the de¬ 


cays) is satisfied [H (IT, 

M^SxloSGevG^) 2 , < 4 > 

where v is the vev that enters the see-saw of Eq. [3] 

The ability of the minimal Type I see-saw model to 
simultaneously explain neutrino masses and the BAU is 
certainly intriguing. However, Vissani observed |l6j that 
the model is incapable of doing so without generating a 
naturalness problem^ Equation [4] is simply incompat¬ 
ible with the conservative naturalness requirement that 
corrections to the electroweak /i 2 parameter of Eq. [l] not 
exceed 1 TeV 2 . With three flavours of hierarchical right- 
handed neutrinos this requires d 

M„. < 3x10’ GeV ■ (5) 

The incompatibility is exemplified in Fig. [T] nowhere at 
v = 246 GeV is it possible to simultaneously fulfil the 
Davidson-Ibarra and Vissani bounds. 

A sensible question is then: in what minimal ways 
can this incompatibility be overcome? Figure |T] suggests 
three conspicuous (but not mutually exclusive) options: 
(1) modify the correction to y 2 , e.g. by restoring su¬ 
persymmetry, or by partly cancelling the correction from 
the heavy fermion loop [HI |T5]; (2) lower the Davidson- 
Ibarra bound, e.g. by considering resonant leptogenesis 
[20] . an alternative mechanism 121] . or by introducing 
new fields which allow an increased CP asymmetry in 
the right-handed neutrino decay; (3) seek an extension 
of the canonical see-saw for neutrino mass, i.e. reduce 
the (possibly effective) v entering the see-saw Eq. [3| 

In this paper we will consider the third option. Specif¬ 
ically we will examine alternative see-saw possibilities 
when the minimal Type I see-saw model is extended by 
a second Higgs doublet <I> 2 . We are motivated by the 
following observation: if the see-saw neutrino mass of 
Eq. [3] is evaluated at v < 30 GeV, then Eqs. [4] and [5] 


1 Naturalness is admittedly an aesthetic requirement of a model, 

and the possibility remains that nature is just fine-tuned. 
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FIG. 1. Bounds on hierarchical thermal leptogenesis as a 
function of v. Shown is the Davidson-Ibarra bound (purple), 
and the Vissani bound (blue). The dashed lines indicate v = 
246, 30 GeV. 


become compatible, as is clear from Fig. |T| Thus we 
expect that two-Higgs-doublet models with right-handed 
neutrinos (^2HDMs) and tan/3 = v\/v 2 > 8, where $2 
is responsible for a tree-level see-saw, can naturally ac¬ 
commodate leptogenesis and neutrino masses. In fact, 
we find that tan /3 > 4 is possible, since the extra scalar 
states can be naturally TeV-scale and the Vissani bound 
can be relaxed. 

Upon examining the J/2HDM scenarios which will suc¬ 
ceed, we rediscover the radiative Ma model [223 as the 
only possibility when V 2 = 0. Otherwise, in models with¬ 
out a significant radiative neutrino mass component, we 
require 0.3 < u 2 /GeV < 60. The potentially small vev is 
made natural by softly breaking a (7(1) or Z 2 symmetry, 
which also automatically results in one SM-like CP-even 
state. One advantage of this scenario is that it can work 
in all 1/2HDM Types, greatly increasing the opportunity 
for model building. 

The paper is organised as follows. In Sec. |TT]we build 
the ^2HDM models of interest, describe the scalar states, 
and briefly review the relevant experimental constraints. 
In Sec. |III| we pay particular attention to naturalness lim¬ 
its on the extra scalars; we verify that a natural Z/2HDM 
of any Type is still allowed by experiment. We discuss 
neutrino masses in Sec. [TV] and leptogenesis in Sec. [V] 
The region of parameter space which naturally achieves 
hierarchical leptogenesis is identified. We conclude in 
Sec. EH 


II. THE v 2HDM MODEL 
A. Lagrangian 


Model 

Ur 

4 

4 

4 

Type I 

$1 



$2 

Type II 

$1 

$2 

$2 


Lepton-specific (LS) 




$2 

Flipped 

$1 

$2 

$1 

$2 


TABLE I. The four models with no tree-level flavour-changing 
neutral currents and allowing for a GeV-scale vev to provide 
the see-saw whilst preserving perturbativity of yt- 


tential (see e.g. Ref. [23]) 

U 2 HDM = rn 2 n + 771.22 - rn \2 ($1$2 + $ 2 $ l) 


($l$i ) 2 + ^($t$ 2 ) 


2 

+ A 3 ($£$ 2 ) +A 4 ($ 2 $ l) 

+ ^f($ t 1 $ 2 ) 2 +($^i) 2 j, (6) 


where all the parameters are real. To explain observa¬ 
tions, at least one of these doublets must obtain a non¬ 
zero vev. We consider m 31 < 0 and a CP-conserving 
vacuum, 


(,f,)o = ( “ 


< '" 2>0 = 71 (»° 


(?) 


where v\ > 0, V 2 > 0, and vf + v 2 = v 2 « (246 GeV) 2 . 

A general 2HDM will have flavour-changing neutral 
currents at tree-level. These can be avoided if right 
handed fermions of a given type (u R , d l R’ e k) couple to 
only one of the doublets JM[ [25]. Although not strictly 
necessary, we will assume that this is realised, and adopt 
the convention that only $1 couples to the u R . In a 
7/2HDM, if we assume this also applies for the v\ ? , then 
there are eight possibilities. As mentioned in the In¬ 
troduction, the see-saw constraint Eq. [3] can be made 
consistent with naturalness and leptogenesis if the vev 
contributing to the see-saw is sufficiently small. Since 
we would like our model to remain perturbative, and al¬ 
ready y t ~ 1 for v « 246 GeV, we anticipate that $2 
obtains the small vev and thus we couple it to the v l R . 
Remaining are four possible ti2HDMs which we refe r to 
by their conventional Types as listed in Table ® The 
Yukawa Lagrangian is then given by 


— C-Y = + Vu<lL^lUR + Vd^L^ldR 

+ Veh^jeR + uJl^ 2^R 
+ ^M N (v R ) c v R + h.c., (8) 


where I, J depend on the model Type, and family indices 
are implied. 


The scalar content of the model contains two doublets 
2 each with liypercliarge +1. For simplicity we con¬ 
sider the softly broken, CP-conserving, ^-symmetric po- 


2 Type I (/2HD\Ts with V2 ~ eV were considered in Refs. |26lf29| . 
We will end up considering i >2 of 0(0.1 10) GeV. 
















B. Scalar masses and mixings 
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Consistency with experiments requires the extra scalar 
states to have masses at least > 80 GeV. In order to 
construct models with potentially TeV-scale scalars with 
a naturally small v 2 , we will consider m 22 > 0 and 
m\ 2 /m\ 2 ^ 1 [26]. This is technically natural, since 
in the limit m\ 2 /m 22 —► 0 a U( 1) or Z 2 symmetry is 
restored if A 5 = 0 or A 5 7 ^ 0 , respectively. 

For m 2 -, < 0, the vevs are given by 


Vl 


V2 



1 


m 12. 


1 + ^A 3 45 m 22 


-Vl, 


(9) 

( 10 ) 


where A 345 = A 3 + A 4 + A 5 . These relations become exact 
when m 2 2 = 0 and when A 2 = 0, respectively. For m 22 
A 345 vf, tan/3 = vi/v 2 ~ m% 2 /m\ 2 . 

There is a useful constraint on ?ti 22 which is derived as 
follows. The minimisation conditions give 


1 


m h + o A 2« 2 = tan /3 m n + -Ai^ 


1 


( 11 ) 


where Ai ^ 2 ~ m ^ ~ (125 GeV ) 2 (see below). In the limit 
to 2 2 X 2 v 2 /2, it can be seen that m 22 is bounded above 
by 


m 22 <^m 2 h tan 2 /3 ( 12 ) 

if mh is to remain negative. Figure[2]illustrates how mh 
deviates from its standard value of —(88 GeV ) 2 as m 22 
approaches this bound. For m 22 above this bound, mh 
very quickly grows to values > v 2 , and v ss 246 GeV is 
only explained by a miraculous balance of mh against 
m 22 / tan 2 /?, which constitutes a fine-tuning. Thus we 
adopt Eq. [12] as a consistency condition. 

The charged scalar and pseudoscalar (neutral scalar) 
mass-squared matrices are diagonalised by a mixing angle 
/3 (a). The neutral mass eigenstates are 


h = pi cos a + p 2 sin a , 

H = p 2 cos a — pi sin a , 

A = r/ 2 sin /3 — pi cos /3, (13) 


where pi = -s/2Re(<E>°) — Vi and rp = v / 2Inr(<l)( ) ). The 
masses are given by 


m 2 = Aiu 2 + O 


m 


12 „,2 


m H = m 22 + + O 


m 


12 


m: 


m A — m 22 + ~ ^s) V l + O 


22 J , 

A 

m 


12 


m 


rrin 


22 I ’ 


m H± = m 22 + 2 X?,Vl 


m 


12 


m 


22 


(14) 


tan/5 
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FIG. 2. Contours of mh/ GeV 2 = —80 2 ,— 70 2 , and so on. 
The shaded region has no solution for mh < 0. 


i.e. the same as in the inert doublet model [23 up to 
corrections proportional to m\ 2 /m \ 2 , which we provide 
in Appendix [A Clearly, if m 22 v 2 , the mass scale of 
extra scalar states is ~ m 22 . 

In the alignment limit cos(a — /3) —>• 0, the couplings 
of h to SM particles become SM-like. We calculate 


cos 2 (a — /3) 


_4 4 

m\ 2 vl 


(Ai — A 345 ) 


4 4 

m 22 m 22 


1 - 


2 mi 


(2Ai — A345 


2 / 2 \ 2 ’ 

( X + 2hfc A 345) 

(15) 


suppressed by the approximate 17(1) or Z 2 symmetry 
(mh/m 22 <C 1 ) as well as the usual decoupling limit sup¬ 
pression ( v\/m 22 <C 1) [30]. Thus the model naturally 
accommodates a SM-like neutral scalar state. 


C. Constraints 

With Mn > 77122 the constraints (and search strate¬ 
gies) for a J/2HDM of given Type are largely identical to 
those for a 2HDM of the same Type, for which there is ex¬ 
tensive literature (see references henceforth). The 2HDM 
potential Eqj 6 ] is subject to a few standard theoretical 
constraints [33]. The necessary and sufficient conditions 
for positivity of the potential in all directions are [311433] 

Al,2 > 0, 

A3 > — \/AlA2, 

A3 + A4 — |A51 > — \/A1A2. (16) 

Vacuum stability of the potential minimum is more diffi¬ 
cult to evaluate. An inequality which ensures a global 
minimum, missing possible metastable vacua, is pre¬ 
sented in Ref. [34] . Tree-level perturbative unitarity of 
scalar-scalar scattering is ensured by bounding the eigen¬ 
values of the scattering matrix Perturbativ- 

ity of the Acan also be demanded [3S]. At the very 
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least these bounds should be implemented at the mass 
scale of the scalar states. In addition they may be de¬ 
manded up to some high scale under the renormalisa¬ 
tion group evolution, which results in non-trivial con¬ 
straints on the parameter space (see e.g. Ref. [55VHTQ . 
Type II, LS, and Flipped 2HDMs are particularly sus¬ 
ceptible to exclusion by such a demand at large tan/I; 
at one-loop, their Yukawa couplings hit a Landau pole 
before Mpi ~ 10 18 GeV when V 2 < 3.6,2.3,3.3 GeV 
(tan/3 > 68,107,75) respectively [42] . These Landau 
poles merely indicate the breakdown of perturbativy. 

The scalar boson discovered at the LHC is to be iden¬ 
tified with the mass eigenstate h. Its couplings have been 
measured to be SM-like, which constrains the ^2HDM to 
lie in the alignment limit cos(a — /3) sa 0, particularly 
at large tan/3, and for Type II and Flipped 2HDMs[^] 
As is evident from Eq. |15[ the alignment limit is auto¬ 
matically preferred in our model due to the approximate 
Z-2 or 17(1) symmetry. Thus we limit the following dis¬ 
cussion on additional experimental limits to those that 
constrain moderate to large tan/3 models very close to 
the alignment limit. 

In Type II and Flipped 2HDMs, the $2 coupling to 
down-type quarks is tan/3 enhanced. The H ± state then 
contributes significantly to radiative B —> X s 7 decay; the 
experimental measurement [491 combined with a recent 
next-to-next-to-leading order SM calculation m bounds 
m H ± > 480 GeV at 95% CL for tan/3 > 2. This bound 
along with the consistency condition Eq. |12| implies V2 /$ 
45 GeV (tan /3 > 5.4) for these z/2HDMs. In the Type 
II 2HDM the $2 coupling to e R is also tan /3 enhanced, 
and the bound on m R ± from B —>• tv decays exceeds the 
radiative bound when tan /3 > 60 (23). 

Direct searches at LEP constrain m R ± > 80 GeV as¬ 
suming decay to SM particles m ■ At the LHC, searches 
for H/A —> tt [52, 53] are particularly constraining in 
the Type II 2HDM. The 95% CL limit rises approxi¬ 
mately linearly from > 300 GeV at tan/3 = 10 to 
m A > 1000 GeV at tan/3 = 60. Such searches can 
also be mildly constraining for the LS 2HDM at mod¬ 
erate tan/3. Searches for if* —> tv jSU 55] cannot 
compete with B —» X s j for Type II/Flipped 2HDMs 
or with H/A -4 tt for the LS 2HDM. However, for 
m H ± < 160 GeV, significant parameter space is ruled 
out in Type I 2HDMs with moderate tan /3. 

The {yv)ijl l L <&2 l ' :l R Yukawa term related to the neutrino 
masses can induce lepton flavour violating decays; these 
are suppressed by the small y v and the right-handed neu¬ 
trino mass scale Mjv > m ,22 ■ The processes of interest 
are l a -4 ly'y, l a -4 3 Ip, and y —> e conversion in nu¬ 
clei (see Ref. [55] for expressions). As well, b —► sl a lp 


3 Note that some of the bounds derived in these papers do not 
apply to the softly broken Zo-syrn metric case, and also do not 
apply to the ^-symmetric case when one of the vevs vanishes 

4 We refer the reader to Refs. [411 |43H48[ for allowed parameter 
space as a function of cos(rv — ft ) and tan ft in all 2HDM Types. 


decays are induced in Type II and Flipped jz2HDMs. In 
practice, lepton flavour violating measurements constrain 
linear combinations of (y u )ij bi- and tri-linears as well as 
the Mjv ; . 

In summary, for moderate to large tan/3 and cos (a — 
/3) ~ 0, experiments are most constraining for the Type 
II and Flipped 2HDMs, with 77122 > 480 GeV necessary 
(implying V 2 < 45 GeV). For Type I and LS 2HDMs, even 
additional scalars with masses down to 80 GeV may still 
have evaded detection. 


III. NATURALNESS 


In the SM, the renormalisation group equation (RGE) 
for the electroweak y parameter (as in Eq. [l]) is domi¬ 
nated by the top quark Yukawa, 


dy 2 _ 1 A 2 2 

d\ny R ~ (47r) 2 Vt ^ 


(17) 


where y R is the renormalisation scale. At low energy 
we measure y 2 ~ —(88 GeV) 2 , and under SM running it 
is apparent that \y\ remains ~ 100 GeV even up to the 
Planck scale Mpi ~ 10 18 GeV. Thus there is no measur¬ 
able naturalness problem in the SM alone; there is no 
fine-tuning of any measurable parameter at a high scale, 
only the cancellation of an unmeasurable bare parameter 
against an unphysical cutoff scale, which should be as¬ 
signed no physical significance. With this understood, it 
is clear that a measurable naturalness problem can only 
arise when dy?/d ln/i# > y 2 . Indeed, this is exactly how 
the Vissani bound in the Type I see-saw model can be in¬ 
terpreted PS HIl. Let us now examine when the Z/2HDM 
encounters such a problem. 

In practice, the naturalness considerations can be di¬ 
vided into two distinct calculations: the influence of TO 22 
on mu, and the influence of M/v on 77122 - These influ¬ 
ences will be considered in turn0 


A. Corrections to mfj 

If m 2 2 7 ? 7 ^ tan 2 /3/2, then m 2 -, sets the mass of the 

observed SM-like Higgs via Eqs. [9] and [14] The one-loop 
RGE for the m 2 X parameter is [28] (see Ref. !H] for a 
recent two-loop calculation) 

^ = jhf [ (4As+ 2A4)W22 + • (18) 

The 0 ( 777 ^) term contains gauge, Ai, and Yukawa con¬ 
tributions, which, as in the SM case, do not induce a 
naturalness problem. However if A 37 are non-zero then a 


5 In the following we ignore the influence of the small y v Yukawas 
on mf 1 , and hence those results also hold in a general 2HDM. 
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naturalness problem is induced for sufficiently large to| 2 ; 
we are interested in when this generically occurs. Even if 
A 3 4 = 0 at some scale, they will quickly be reintroduced 
by gauge interactions at one-loop. Their one-loop RGEs 
are given by 


d\ 3 
din Hr 
d\± 
din hr 


1 

( 4 - 7 t) 2 

1 

( 47 t) 2 


4 {dr — 2 gyg\ + 3 g 2 ) + ... 
[Z 9 W 2 + 


(19) 


where g 2 {mz) ~ 0.43 and gy(mz) « 0.13 are the gauge 
couplings and the ellipses contain terms multiplicative in 
A 3 , 4 , terms proportional to A 2 , and terms related to the 
Yukawas. Let us ignore those effects for now and return 
to them later. Note that ignoring the contribution from 
A 2 is equivalent to assuming A 5 < 0.2, so that its contri¬ 
bution is subdominant to the gauge couplings. Typically, 
one would expect 

IW(l« r )I ^ (4^)2 4 ~ ^ 9 y 92 + 3^2) i 

\^4,(Hr)\ > ( 4^)2 ^9y92i (20) 

and thus 

> ^4 0* 9 y + 0^2) m 22- (21) 


dm{i 

d In h r 


This lower bound is of the same order as the two-loop 
pure gauge contribution Ell- 

Equation 21 represents a conservative bound on the 
running of the mf 1 parameter above the scale ~ m 22 - 
Naturalness demands that this running not be signifi¬ 
cantly larger than the value measured at a low scale, 
|mn| « 88 GeV. A very conservative naturalness bound 
is therefore 


7^4 ( 3 9y + 9s 2 4 ) m \ 2 < 1 TeV 2 , (22) 

=> m 22 < 1 x 10 5 GeV. (23) 


Alternatively, we can try to bound a quantity which 
measures the fine-tuning in to)-, at some high scale A/,. 
A typical quantity is wnm 


m n(A/i) drn^j 0) 
m^(0) dml^An) 


(24) 


which compares percentage variations of two (in princi¬ 
ple) measurable parameters. Let us now estimate how 
such a bound might constrain 77122 - 

For simplicity, and anticipating that the 77122 scale 
is not far above the electroweak scale, we will evolve 
the dimensionless parameters using the (i/)2HDM RGEs 
from the mz scale. First, the one-loop gauge coupling 
RGEs [23, 122) can be solved analytically. Upon sub¬ 
stitution into the A 3,4 RGEs (Eqs. 19), and considering 


only the pure gauge contribution, the A 3,4 running can 
be solved for given initial conditions. For simplicity we 
take \ 3 {m 22 ) = X^Tn. 22 ) = A 3 , 4 ( 77122 ) and consider it a 


free parameter. Next we solve Eq. 18 for 771)4 (hr) with 


the initial condition 77744(77722) = —(88 GeV ) 2 (neglecting 
any RGE evolution of to 22 ). With these simplifications 
977744 (0)/<9 t 77 4 4 (A ft ) = 1 and the fine-tuning measure is 
given simply by A (Aft) = 1 777.^4(A*,)/(88 GeV) 2 |. 

Note that in setting the initial condition 777^4(777-22) = 
— (88 GeV ) 2 we have implicitly assumed that m | 2 -C 
to); tan 2 / 3/2 (see Eq. 12 and Fig. HI). This is conser¬ 
vative for negative 777^1 , since 1 777)4 (77722)! shrinks as 
to 22 / tan 2 /3 —> to);/ 2 and the naturalness constraint 
would become more stringent. In some circumstances 
we will obtain naturalness bounds on m 22 which ex¬ 
ceed to 2 tan 2 / 3 / 2 , which just indicates that the natural¬ 
ness constraint is weaker than the consistency condition 
Eq.[l 2 } 

In Fig. [ 3 ] we show A = 10 and A = 100 contours as a 
function of Aft and A 3 , 4(77722)- These represent natural¬ 
ness upper bounds on 77122- The cusp-like structures of 
apparently low fine-tuning in 77744 occur when 77744 runs 
negative before turning and passing through 777)4 = 0 , 
which just corresponds to a fine-tuning in (A 3 , 4 , Aft). A 
stringent naturalness constraint is obtained by demand¬ 
ing A < 10 at Aft = Mpi ; from Fig. [ 3 j it is clear that this 
implies 


77122 few x 10 3 GeV. 


(25) 


If any new physics comes in below Mpi then the running 
of TO 44 could change, and these bounds do not apply. If 
that is the case then it is more appropriate to consider 
Aft at the scale of the new physics, which weakens the 
bound, as is clear from Fig. [3j In the 772 HDM this new 
physics scale is the right-handed neutrino scale Mat, after 
which the right-handed neutrinos can contribute to the 
running of 77744 through to 22 at one-loop. 

We have so far ignored the RGE contributions from 
possibly large Yukawas. There are two situations in 
which the Yukawas play a significant role. The first is 
in Type II, LS, and Flipped t/2HDMs with tan/3 large 
enou gh su ch that an early Landau pole is induced (see 
Sec. IIC), and the second is in Type II and Flipped 


z72HDMs with moderate to large tan/3 when the pure 
Yukawa term ± 12 t/ 2 t// induced by a quark box dia¬ 

gram contributes significantly to the A 3,4 RGEs (Eqs. 191. 
In Fig. [4] we show how the A = 10 contours change as a 
function of V 2 in an example Type II I/2HDM. For this 
Figure we have numerically solved the full set of one- 
loop RGEs [JB including the top/bottom/tau Yukawas, 
taking the following values at the scale Mz : <? 2 = 1.48, 
Ai = A 2 = 0.26, and y t = 0.96/sin/3, yb = 0.017/cos/3, 
y T = 0.010/cos/3. Comparing to Fig. |3j it can be seen 
that the pure Yukawa term has a noticeable effect when 
V 2 /$ 20 GeV. It is also apparent from Fig. [4]that nearing 
V 2 ~ 3.6 GeV (below which a Landau pole is induced be¬ 
fore Mpi) can act to degrade or improve the naturalness 
bound. The V 2 = 3 GeV bound in Fig.[4]shows the effect 
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V,4( m 22) 

FIG. 3. Contours of the fine-tuning measure A (A/,) = 
10 (100) in black (grey) as a function of (top) for 
A 3 , 4 ( 77 x 22 ) = 0.0, 0.1,—0.01 (solid, dashed, dotted), and (bot¬ 
tom) A 3 , 4 ( 0122 ) for A h = 10 18 ,10 12 ,10' GeV (solid, dashed, 
dotted). See the text for the assumptions that accompany 
this plot. 


of hitting the Landau pole at ~ 10 9 GeV. We note that 
this only signals the breakdown of perturbation theory, 
and of our one-loop RGEs; we cannot calculate 77122 (^ 1 ?) 
above this scale though it is perfectly possible that the 
theory remains natural. 


In a repeated full one-loop RGE analysis we found that 
the Flipped v 2HDM gave essentially the same results as 
the Type II n2HDM in Fig. [4j and there was no notice¬ 
able Yukawa effect in the LS ^2HDM until the Landau 
pole was reached. Thus we find that the stringent natu¬ 
ralness bounds of Eq. 25 and Fig. [3] are applicable at all 
times in the Type I Z/2HDM, for u 2 > 2 GeV in the LS 
J72HDM, and for v 2 > 20 GeV in the Type II and Flipped 
z72HDMs. Otherwise Yukawa effects must be taken into 
account. Either way, the important point is now clear: 
a TeV-scale to 2 2 can be both completely natural and, as 
was discussed in the previous subsection, is experimen¬ 
tally allowed in all Z/2HDM Types. 




b,4( m 22) 


FIG. 4. Contours of the fine-tuning measure A(A h) = 10 in 
an illustrative Type II i/2HDM (see text) for different values 
of i> 2 . The solid lines for V 2 = 20 GeV match approximately 
onto the solid lines in Fig. [3] 


B. Corrections to m| 2 

Let us now consider the influence of the right-handed 
neutrinos. The one-loop RGE for m| 2 is HUES 

^ = (iy[-4TV[^M + 0 «)]. (26) 

A conservative naturalness bound is obtained by bound¬ 
ing the running as we did in Eq. |23[ 

^^Tr [y v 'D 2 M y} / \ < A 2 ound , (27) 

where taking Abound. = 1 TeV gives the Vissani bound 
on M Nl of Eq. [5] [T7] . However, now we are bounding 
corrections to rn% 2 rather than , which may be TeV- 
scale. Thus, depending on the mass of the extra scalars, 
it is possible that we can sensibly take A bound. > 1 TeV, 
in which case the the naturalness bound is somewhat re¬ 
laxed; in Fig. [5] we show a Relaxed Vissani bound for 

Abound = min(10 TeV, 10 \jnn\ tan 2 /3/2), where we have 
kept in mind the consistency condition Eq. |12| The Vis¬ 
sani bound still represents the unnatural area of param¬ 
eter space if TO 22 is closer to 100 GeV. 
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FIG. 5. Bounds on the 1 / 2 HDM as a function of V 2 . Shown 
(as labelled) are the Davidson-Ibarra bound, the Vissani and 
Relaxed Vissani naturalness bounds, and the areas of param¬ 
eter space with strong A L = 2 scattering washout. The Type 
II and Flipped t/2HDMs are excluded by B —> X s 7 for val¬ 
ues of V 2 greater than indicated by the grey dashed line (see 
Sec. 


the 


II C I . The grey dotted lines indicate the V 2 below which 


Yukawas hit a Landau pole before Mn 1 in the Type II, 


Flipped, and LS t/2HDMs riglit-to-left. 


As before, we could instead bound a quantity which 
measures the fine-tuning in m| 2 at some high scale A^. 
In this case, the fine-tuning measure of Eq. [24] is 


A(A h) — 1 + 


1 T,i < j(yv)ij M Hyt)ji ln OV A h) 


47T 2 


m 


22 


(Mj) 


(28) 


Taking 7712 2 (Mj ) ~ 1 TeV and demanding A (Mpi) < 10 
gives a similar bound to Vissani (Eq. 27 with A b oun d = 


1 TeV). Note that there is no naturalness bound on Mn 
in the t/„ —► 0 limit. This is the technically natural 
limit corresponding to an enhanced Poincare symmetry 
in which vr decouples from the theory Ell¬ 
in summary, there are up to three scales in the 
J/2HDM: v, m 22 , and Mjy We have described the con¬ 
ditions under which v 2 (or m\{) is protected from m 22 , 
and m\ 2 from Under such conditions it follows that 
to 2 1 is also protected from and the model is entirely 
natural. 



10 y 

10 s 
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10 6 

10 5 
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10 3 


FIG. 6. As in Fig. [l] but for the Ma model. The Vis¬ 
sani and Relaxed Vissani bounds are evaluated at m 2 2 = 
100,1000 GeV respectively. The Davidson-Ibarra bound and 
strong A L = 2 washout region are shown for 77122 = 500 GeV, 
though they are only mildly sensitive to 77722- 


A. ml 2 >0, A 5 = 0 


If to 2 2 > 0 then the 17(1) lepton number symmetry is 
softly broken, i.e. the breaking does not force us to insert 
a non-zero A 5 term in order to introduce a divergent coun¬ 
terterm, and it is consistent to consider ra^ > 0 , A 5 = 0 . 

In this situation the neutrino mass matrix is given by 
the see-saw formula, 


™ _ 2 ^-,-1 T 
m v — 2 Vv 


1 v“ 
tan 2 (3 2 


yuT> M y v > 


(29) 


TO 22/ TO 12 for m 22 


> A 345 U 2 (see Eq. 10 1 . 


where tan /? 

The analogous Davidson-Ibarra and Vissani bounds 
are given by the standard Eqs. [4] and [5] with the replace¬ 
ment v —> v 2 . These bounds are depicted in Fig. [5j If 
V 2 < 30 GeV (tan/3 > 8) then both bounds are satisfied. 
As well, as discussed in Sec. IIIIl if 777-22 is TeV-scale the 


Vissani bound can be relaxed and the required CP asym¬ 
metry needed to reproduce the BAU via leptogenesis may 
be naturally achieved for u 2 < 60 GeV (tan/3 > 4). In 
the Type I 7/2HDM, v 2 can be naturally <C GeV. Other¬ 
wise, requiring a perturbative theory up to M^ 1 restricts 
v 2 > 1 GeV [v 2 > 2 GeV) for the LS (Type II/Flipped) 
722HDM in the parameter space region of interest, as de¬ 
picted in Fig. [5] 


IV. NEUTRINO MASSES 


If 


to 2 2 = 0 and A 5 = 0 then a 17(1) lepton number 
symmetry can be defined and neutrinos remain massless. 
Let us now consider turning each non-zero in turn. 


B. 7775*2 = 0, A 5 yf 0 

In this situation u 2 = 0 and a Z 2 symmetry remains 
unbroken; this is the scenario of Ma 22j. The model 
yields a radiative neutrino mass and a dark matter can¬ 
didate. This is only possible in the Type 1772HDM, since 













in any other Type the unbroken Z 2 forbids a Dirac mass 
term for any charged fermion coupling to $ 2 - Note that 
the limit A 5 —> 0 is technically natural, since in that limit 
the 17(1) lepton number symmetry is reinstated. 

If My m% 2 ,v 2 the radiatively induced neutrino 
mass matrix is 


:, the rates 


rn = ~(yly v )uMi, 

(32) 

17 T 2 

(33) 

H « 

Mpi 




V )kj W 

~2 W k 8^2 


In 


2 Ml 


L H 



are typically rescaled and expressed in terms of an effec¬ 
tive neutrino mass toi and an equilibrium neutrino mass 
m*, 


The analogous Davidson-Ibarra and Vissani bounds are 
given by the standard Eqs. [4] and [5] with the intuitive 
replacement v 2 —> (ln"[ 2 / (mjj + m^)] — l). 

This assumes that there is no fine-tuning in the complex 
Uv parameters to reproduce the observed neutrino masses 
(see Appendix [B] for details). These bounds are depicted 
in Fig. [ 6 j where the Davidson-Ibarra bound has been 
evaluated at TO 22 = 500 GeV as an illustrative exam¬ 
ple (the bound is only mildly sensitive to 17122 )- We find 
that the Ma model with A 5 < 0.5 can naturally achieve 
the required CP asymmetry to reproduce the BAU via 
hierarchical leptogenesis)^] 


C. m\ 2 >0, A 5 7 ^ 0 

In this case both the tree-level see-saw and the ra¬ 
diative mechanism will contribute to the neutrino mass. 
Both contributions are calculable, and either might dom¬ 
inate. Note that it is still technically natural to take 
A 5 —> 0 in this case, since it restores a softly broken U( 1) 
symmetry. In other words, the A 5 RGEs to all orders 
will be multiplicative in A 5 , indicative of the fact that the 
soft-breaking term can only generate finite U (l)-breaking 
corrections. 


V. LEPTOGENESIS 


The observed BAU is achieved analogously to stan¬ 
dard hierarchical thermal leptogenesis US; the out-of¬ 
equilibrium CP-violating decays of the lightest right- 
handed neutrino N\ —► 2 create a lepton asymmetry 
which is transferred to the baryons by the electroweak 
sphalerons above T ~ 100 GeV. 

The details of the leptogenesis are largely defined by 
the decay parameter 


K = 


Tp 

H\t=m 1 


TO 1 
TO* ’ 


(31) 


comparing the rate for decays and inverse decays to the 
expansion rate at the time of departure from thermal 


6 A similar observation was made in a recent paper m- 


TO 1 

TO* 


(ylyv) iiv 2 

2 Mi 


1.1 x 10 “ 3 eV 


( ___V 

V246 GeV/ ’ 


(34) 


where v is the vev that enters the see-saw Eq. [ 2 ] In 
the Z/2HDM with A 5 = 0 (with m\ 2 = 0, A 5 7 ^ 0) the 
analogous definitions make the replacement v 2 —► v 2 
(v 2 —> v 2 ^ (In [ 2 M 2 N J(m 2 H + m\)\ — l)). Note that 
for the scenarios we are interested in (e.g. v 2 <C v), to* 
is smaller than its usual value in standard leptogenesis. 

When only decays and inverse decays are considered, 
leptogenesis for given K proceeds exactly as in standard 
hierarchical thermal leptogenesis (see e.g. Ref. [53] for a 
review). In the weak washout regime K <C 1, the baryon 
asymmetry strongly depends on the initial asymmetry 
and the initial N\ abundance, with N\ decays occuring 
at T <C Mi. The strong washout regime K 1 is inde¬ 
pendent of the initial conditions, and the asymmetry is 
generated as the Ni fall out of thermal equilibrium. 

The 2f>2 scatterings with A L = 1 (see e.g. Ref. [53] ) 
provide a correction to the simple decays plus inverse 
decays picture; they act to increase N± production at 
T > Mi and contribute to washout at T < Mi. In 
standard hierarchical thermal leptogenesis, the scatter¬ 
ing contributions involving the top quark and the gauge 
bosons are roughly equal. In the present model the gauge 
boson contribution is the same as in the standard sce¬ 
nario. However, by construction, the $2 involved here in 
leptogenesis does not couple directly to the top quark, 
and thus the usual s-channel (Nl •<-> tq) and t-channel 
(Nt O Iq, Nq -f-)- It) scattering contributions do not oc¬ 
cur. Instead, at large tan /3 they can be replaced by 
the analogous contribution from other charged fermions, 
i.e. the bottom quark in Type II and Flipped i^HDMs 
and/or the tau lepton in Type II and LS z^2HDMs. A 
large tau lepton Yukawa will also introduce new s-channel 
(V4 >2 O r$ 2 ) and t-channel (iV$ 2 ^ ‘lud^) 

scattering contributions. All of these processes are pro¬ 
portional to {yly v )n and hence Mifhi/v 2 , with the ap¬ 
propriate ^2HDM replacement for v 2 . Therefore, they 
scale with the decays and inverse decays so that they 
represent only a minor (but obviously important) depar¬ 
ture from the standard leptogenesis scenario. 

The 2 o 2 scatterings with A L = 2 mediated by 
the right-handed neutrinos (4> 2 Z •H- $ 2 ^ < f> 2$2 -H- ll) oc¬ 
cur as they do in the standard scenario. These pro¬ 
cesses are proportional to Tr[(j/ 1 / j/J)(y !y y/’)l] and hence 
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Mini 2 In' 1 where m 2 = Yh m i is the neutrino mass scale 
> 0.05 eV. Comparing this rate to the decay/scattering 
rates oc Mifhi/v 2 , it is easy to see that after mak¬ 
ing the appropriate Z/2HDM replacement for v 2 , e.g. 
v 2 —vl <C v 2 , these scatterings will become compar¬ 
atively more important than in the standard case. For 
T < Mi /3 the thermally averaged A L = 2 scattering 
rate is well approximated by [(JS] 


T AL=2 _ T /246 GeV\ 4 / m \ 2 

H ~ 2.2 x 10 13 GeV V v ) ^0.05 eVJ ’ 

(35) 


where the previously described ^2HDM replacements for 
v 2 hold (see Appendix | b|). In Fig. [ 5 ] (Fig. [6|, we show the 
region in the A 5 = 0 (m\ 2 = 0, A 5 7 ^ 0) Z/2HDM where 
these scatterings are still in equilibrium at T < Mi/3^ 
This is the region where strong A L = 2 scatterings 
can potentially wash out the generated asymmetry, de¬ 
pending on the details of the leptogenesis (e.g. in a 
weak washout scenario with Ni decays at T <C Mi this 
washout may be avoided). Demanding that the scatter¬ 
ings fall out of equilibrium before sphaleron freeze-out 
at T ~ 100 GeV provides a lower bound v 2 > 0.3 or 
A 5 > 10 5 ; this is represented by the strong A L = 2 
scattering washout regions in Figs. [5] and [6] respectively. 
We note that this calculation has been performed in the 
context of a perturbative theory. This is reliable for 
the Type I Z/2HDM but not for Type II, LS, or Flipped 
^2HDMs with sufficiently small v 2 , when perturbativity 
breaks down. 

Putting this all together, we can now read off from 
Figs. [5] and [ 6 ] the regions of parameter space which can 
achieve natural hierarchical thermal leptogenesis. For 
i/2HDMs with m\ 2 > 0 and A 5 = 0, we find 10 3 GeV < 
M n 1 < few x 10 7 GeV is viable for Type I i/2HDMs, and 
10 4 GeV < < few x 10 7 GeV for all other Types if 

they are to remain perturbative. For the Ma model with 
m 2 2 = 0 and A 5 7 ^ 0 , we find viable parameter space for 
10 3 GeV < M Nl < 10 s GeV and 10 " 5 < A 5 < 0.5 

Lastly we note that the lightest scalar state in the Ma 
model is stable)^] It is therefore possible that this state, if 
it is neutral, constitutes some or all of the observed dark 
matter. During the leptogenesis epoch, >I> 2 is produced in 
abundance in Ni decays. Overproduction of dark matter 
is of no concern as long as $ 2 efficiently thermalises at 
or below the temperatures when Ni decays occur, which 
suggests m 22 <C M Nx . In this case the lightest state is a 
thermal relic dark matter candidate. 


VI. CONCLUSION 

The minimal Type I see-saw model is unable to explain 
neutrino masses and the BAU via hierarchical thermal 
leptogenesis without ceding naturalness. The main con¬ 
clusion of this paper is the observation that a second 
Higgs doublet can avoid this problem. These ^2HDM 
models provide a natural solution by reducing the (pos¬ 
sibly effective) vev entering the see-saw formula. This 
can be done radiatively, or by having the second Higgs 
doublet provide a tree-level see-saw with a small vev v 2 , 
kept natural by softly breaking a U( 1) or Z 2 symmetry. 

The models naturally accommodate a SM-like Higgs 
and predict the existence of approximately TeV-scale ex¬ 
tra scalar states in order to remain natural. We redis¬ 
covered the radiative Ma model as the only possibility 
when v 2 = 0; in that case we found 10 3 GeV < M.y, < 
10 8 GeV and 10 -5 < A 5 < 0.5 could simultaneously ex¬ 
plain neutrino masses and the BAU via leptogenesis while 
remaining natural. The v 2 > 0 models require tan f3 > 4; 
we found 10 3 GeV < M.y, < fewx 10 7 GeV was viable for 
Type I z/2HDMs, and 10 4 GeV < M Nl < few x 10 7 GeV 
for all other Types if they are to remain perturbative. 
The interesting areas of parameter space are well sum¬ 
marised in Figs. [5] and [ 6 ] 
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7 A similar plot to Fig. [H] appears in Ref. Eg in the context of 
the Type I z^2HDM with i >2 > 0. We are not aware of any plot 
similar to Fig. [6] in the literature, though see Refs. 16614701 for 
leptogenesis studies at points in the Ma model parameter space. 

8 The lifetimes of the heavier scalar states are governed by mass 
splittings A via T ~ G^A 5 /(10 2 7r 3 ). In the parameter space of 
interest, one can check that A is typically already large enough at 
tree-level so that lifetimes remain well below 0(1 s) and therefore 
do not disturb big bang nucleosynthesis. 
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Appendix A: Squared scalar masses at 


To order mf 2 / m 225 the scalar masses Eqs. 14 are given by 


Ai + 


‘12 


2A345 — Ai — 2^77 A1A345 




m\ 


m 2 H ± 


( X + 2^ (A345 - 2 Ai)) (l + 2 ^A 3 45 ) 


1 + A; 


tf, 1 _ 2^r (2A345 - 3 A 2 ) + (a§ 45 + 3 A 2 a 34 5 - 6AiA 2 ) 


345 


2?7l22 


m 22 


‘22 


1 + (A345 _ 2A5) 


‘22 


2 m | 2 


i 1 + 2^ ( A 345 - 2Al)) (l + 2^A 3 45 


+ 


mf, 1 + 2m:l : ( A 345 + A 2 - 2A 5 ) 


2 m^ 


-A 345 


1 + Aj 


A 1 + 


‘12 


(A2 + A3) 


2?Tl22 


1 + 


-A345 


(AI) 


Appendix B: Bounds for more general m„ 


We consider the see-saw Lagrangian as in Eq. [2] and a 
neutrino mass matrix of the form 


m. v = l —y u V M 1 V f{M )y ) 


T 
v 5 


(Bl) 


where V x = diag(xi, x%, £ 3 ). Note that in the ^2HDM 
with A 5 = 0 we have f{Mj) = v^/v 2 , and in the Ma 
model with Mjv m 22 we have 




Mf 


- 1 


{m 2 H +m 2 A )/2 
Following Casas-Ibarra eh, it is possible to write 

yj 2 + I I _ I 

?/ — __ TUT) 2 RT) 2 7) 2 

yv — u M ^ f(M ) 5 


(B2) 


(B3) 


where R is a (possibly complex) orthogonal (RR T = 
f? T f? = I) matrix. The Vissani bound on each right- 
handed neutrino mass becomes mi 


1 2 


M .1 


An 2 v 2 f(Mj) S m< l i2 «| 2 < 1 TeV ' 
=> Mjvj < 3 x 10 7 GeV x f{M N Js. 


(B4) 


r 


The CP asymmetry for hierarchical neutrinos [14] be¬ 
comes 


l e i I = 

< 


6 M 1 lui[(R^V m (RD f ^ M) R T )V m R*)i 1 ] 


87 T V 2 

6 m 3 M Nl 


(RAT> m R) 11 


87 t v 2 min[/(Mj)] 


(B5) 


where the approximate inequality holds for max(|f?.jj|) < 
1. For larger max(| R t j |) the inequality can be exceeded, 
but this corresponds to a fine-tuning (see Ref. [T7]). With 
this caveat, the Davidson-Ibarra bound for min[/(Mj)] = 
f(M Nl ) therefore becomes 

M Nl > 5 x 10 s GeV x f(M Nl ). (B 6 ) 

The A L = 2 scatterings are proportional to puna 

1 


J2Re[(ytyv)ij(yty> 




hJ 


Mi Ma 


= Tr [{y v V^yl){y u V^yl)^] 

= M [D m {RV~l M) R T )V m {RV~l M) R T ) t] 


< 1 - 


- 0 

m 


v A min[/(Mj )] 2 ’ 


(B7) 


where the last line is an exact equality for /(Mi) = 
/(M 2 ) = /(M 3 ), an exact inequality when R is real, and 
an approximate inequality (as indicated) if R is com¬ 
plex with max(|f?y |) < 1. Again, for larger max(|i?jj|) 
the inequality can be exceeded. Then for min [f(Mi)\ = 
/(Mjvj), the A L = 2 scattering Eq. 35 becomes 


Tal=2 

H 


< 


T 


1 


2.2 x 10 13 GeV f{M Nl ) 2 


m 

0.05 eV 


■ (B 8 ) 
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